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Abstract. In a recent experiment, Hirjibehedin et al. [Science 317, 1199 (2007)] 
performed inelastic tunnelling spectroscopy of a single iron atom absorbed on a 
nonmagnetic substrate. The observed steps in the differential conductance marked 
the spin excitation energies. In this paper, we explain observed nonmonotonicities in 
the differential conductance by a nonequilibrium population of the atom spin states. 
Furthermore, we predict super-Poissonian current noise due to this nonequilibrium 
situation. We argue that the remarkable absence of nonequilibrium features at certain 
conductance steps indicates the presence of an anisotropic relaxation channel. 
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1. Introduction 

Inelastic scattering processes provide a convenient tool to study the excitation spectra 
of various systems. By using inelastic electron tunnelling spectroscopy, one can access 
vibrational excitations of ensembles of molecules in metallic tunnel junctions [1, 2] or 
single molecules in scanning tunnelling microscope (STM) geometries [3]. The excitation 
energies reveal themselves as steps in the differential conductance whenever a new 
inelastic transport channel opens up. For the explanation of the measured signals an 
equihbrium distribution of the molecule states was implicitly assumed. Further studies 
of molecular vibrations were performed using H2 molecules [4], Cqq molecules [5] in 
mechanical break junctions and suspended carbon nanotubes [6, 7, 8]. 

Recently, the investigation of magnetic properties and interactions on an atomic 
level became possible due to the advent of spin inelastic electron spectroscopy [9, 10, 
11, 12, 13, 14]. Here, single magnetic atoms absorbed on a nonmagnetic substrate 
were contacted using an STM tip. Describing the atom in terms of a localized spin, 
Hirjibehedin et al. [12] related the positions of the conductance steps to the energy 
associated with transitions between different eigenstates, while the relative step heights 
depend on the matrix elements of the spin operator. A more complete theoretical 
description based on perturbation theory in the tunnel coupling [15, 16, 17, 18, 19] still 
relies on the assumption of equilibrium occupations. 

While the above studies could explain the conductance steps assuming the atom spin 
to be in thermal equilibrium with the substrate, nonmonotonic features clearly present 
in the experimental results of [12] were not addressed. Conductance overshoots due to 
nonequilibrium occupations together with their relaxation by spin-phonon interactions 
have been discussed in [20] for cotunnelling through a quantum dot. A similar behaviour 
was found in [21, 22, 23] who additionally studied the low-temperature nonequilibrium 
logarithmic Kondo enhancement of this overshoot. While the Kondo effect is relevant 
for transport through a single Co atom studied in [13, 14], it is not important here. 
Nonequilibrium effects have been considered in [19] for spin-transfer torque on a single 
atom coupled to ferromagnetic substrates and tips. In [24], the nonequilibrium current 
and current noise through a single molecular magnet was analyzed in the charge 
fluctuation regime. In this paper, we explain the experimental results by calculating the 
nonequilibrium occupations together with a spin- dependent relaxation channel using a 
master-equation approach. We, furthermore, predict an enhanced Fano factor indicating 
super-Poissonian current noise as a clear sign of a nonequilibrium situation. 

2. Model 

We model the experimental setup of [12] as two reservoirs of noninteracting electrons 
coupled by a tunnel barrier with an embedded spin. Hence, the Hamiltonian describing 
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Table 1. Eigenenergies E„i and eigenstates |m) of the spin Hamiltonian (2) in the 
basis \ Sz)z of the eigenstates for a magnetic field apphcd in the z and the x direction, 
respectively. 
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the system is given by 

H = ^ Hr + i^spin + -f^tun- (1) 

r 

Here Hj. = ^rkO^lka^rka Hiodels the two electrodes as reservoirs of noninteracting 
electrons with constant density of states pr and electrochemical potential fir- The 
operator a^^^^ creates an electron in lead r = L, R with momentum k and spin a. The 
local spin is described by 

Fspin = -DS', + E{Sl - si) + gp^B ■ S, (2) 

where the z axis is the magnetic easy axis of the atom in its coordination environment. 
For an 5" = 2 iron atom on CU2N, the best fit to the experimental results in [12] gives 
an uniaxial anisotropy D = 1.55 meV, a transverse anisotropy E = 0.31 meV and a 
(jf-factor of (? = 2.11. In table 1, we summarize the eigenenergies and eigenstates of the 
spin Hamiltonian for two different choices of the external magnetic field. Finally, the 
tunnelling Hamiltonian is given by the Appelbaum Hamiltonian [25] 

-f^tun = jrr'a'l]f^cr~~~^ C^r'kV', (3) 

rr'kk'fTcr' 

with cr denoting the Pauli matrices, which describes an exchange interaction between 
the spin of the tunnelling electron and the local spin. We neglect direct tunnelling 
through the barrier not involving the localized spin as it only gives rise to a bias- 
independent elastic background to the differential conductance. Interference terms 
between direct and exchange tunnelling do not appear in the total current and shot 
noise since the contributions from spin-up and -down electrons cancel out each other 
for nonmagnetic electrodes. The above model has been studied extensively to describe 
molecular magnets. The Kondo effect induced by the transverse anisotropy [26, 27, 28] 
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as well as Berry phase effects [29, 30] and the current-induced switching of the molecule 
spin [31] has been discussed. 

We parametrize the couplings jVr' through the sum J = jLL+jRR and the asymmetry 
a = (jLL - jRR)/(iLL + Jrr), 1-6., j^l = (1 + a)VV4, jIr = (1 - a)VV4 and 
Jlr = Jrl — (1 — a^)</^/4. While the couplings jlr and jrl are responsible for the current 
through the atom, which may be accompanied with a spin excitation or disexcitation, 
the couplings jll and jrr do not contribute to the current but give rise to a transport- 
induced relaxation mechanism for the local spin only. 

The dynamics of the system is governed by a generalized master equation for the 
probabilities Pm to find the spin in one of its eigenstates |m) with energy Em, 

^(t) = J2 / dt'Lmm'it - t')Pm'{t') . (4) 

In the stationary limit, Pm = Pm{t) is independent of time and only the time- 
integrated kernel Lmm' = J^^dt'Lmm'i—t') is needed. Its matrix elements are Lmm' = 
Wmm' — ^mm'Wm- Here, Wmm' are the Fermi's golden rule transition rates, 

Wmm' = ^ 27r \jrr' \ "^prPr' \{m\SJ\m')\^C, {fir " Pr' " \am' ) , (5) 
rr'a 

where ({x) = — e"^/^'^'^"^^), and A^m' = — Em'- The sum runs over the lead 
indices r, r' = L, R and the spin directions a = x,y, z. The elements Wm follow from 
^m ^ram' = which guarantees the conservation of probability. 

To compute the current and current noise, we employ the formalism of full-counting 
statistics adopted to system that can be described by rate equations [32, 33]. To this end, 
we introduce the matrix W^^, which is obtained from Wmm' by multiplying each term 
in the sum (5) with a factor e*-^ if r = L, r' = R, e~*^ if r = R, r' = L, and otherwise, 
where x is called a counting field. Furthermore, we define L^^, = W^^, — Smm'Wm 
(note that Wm does not contain the counting field x)- The smallest eigenvalue of L^^, 
defines the cumulant generating function S{x), from which we can obtain the average 
current / and the current noise S by performing derivatives with respect to the counting 

field, I = -ie {dS{x)/dx\=o and S = (-ie)^ {d'S{x)/dx%=o- 

Although the full-counting statistics formalism to compute the current and noise 
is very compact and elegant for the calculation, we introduce, in addition, an 
equivalent formulation for the average current, that offers a more transparent basis 
for distinguishing equilibrium from non-equilibrium effect. It is easy to show that the 
average current can also be written as 

'^"'^ Pm'. (6) 



W^ , 

, mm' 



The derivatives i(dW^„//dx)lx=o are the current rates. Non-equilibrium effects of the 
current- volt age characteristics enter via the non-equilibrium probability distribution Pm, 
that is obtained by solving the master equation, Eq. (4). These non-equilibrium effects 
would be neglected if one replaced the Pm by an equilibrium probability distribution, 
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Figure 1. Differential conductance in units of Go as a function of bias voltage for 
different values of the asymmetry parameters a. Arrows indicate the position of 
the excitation energies. Parameters are = 7T and T = 0.5 K corresponding to 
the experimental values of [12]. The corresponding eigenenergies and eigenstates are 
summarized in table 1. 



= exp{—Em/kBT)/J2m'^^'Pi~^m'/kBT), i.e., for low temperature Pq = 1 for the 
ground state and Pm = for the excited states m ^ 0. 

3. Results 

In the following we discuss the influence of a nonequilibrium spin occupation on the 
transport properties for the system parameters of the experiment [12]. In figure 1, we 
show the differential conductance in the presence of a strong magnetic field = 7T 
along the easy axis for different values of the asymmetry parameter a in the absence of 
the phenomenological relaxation (12), see discussion below. For very large asymmetries, 
a — )■ 1, there are fiat plateaus between the conductance steps. In this limit, the coupling 
constant jin for processes that drive the atom state population out of equilibrium is 
much smaller than Jll for processes that let the system relax to thermal equilibrium 
with the left electrode. Therefore, as in [1], nonequilibrium effects are absent, and 
the resulting conductance curve is identical to the one obtained in [15]. For smaller 
asymmetries, the situation is different. The height of the conductance steps at the 
excitation thresholds is increased. Beyond the threshold voltages, the differential 
conductance shows a slow powerlaw decay towards its value for the equilibrated system 
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again. This overshooting behavior is observed for most of the steps in the experiments 
of [12]. While the couphng of the adatom to the substrate is fixed in experiment, 
the couphng to the STM tip can be controUed by changing the tip-atom distance. This 
corresponds to changing the total coupling J and thereby the total tunnel current as well 
as the asymmetry a and thereby the nonequilibrium effects. In a recent experiment using 
a magnetic STM tip [34], it was confirmed that by decreasing the tip-atom distance and 
therefore increasing the current through the system, the nonequilibrium effects became 
more pronounced. 

Before we discuss this for the system at hand, we illustrate the mechanism that 
leads to this conductance behaviour explicitly for the simpler model of a local spin-1/2 
with Zeeman energy B, symmetrically coupled to the electrodes, at zero temperature. 
Transport takes place by either spin-flip or spin-conserving transitions. The latter 
contribute to the current as 

^sc = TreliLRppLPReV, (7) 

independent of the probabilities and to find the spin in state up and down, 
respectively. The differential conductance, measured in units of Gq = 47ie'^S{S + 
l)bLRpPLPR is Gsc = Go/3. Therefore, nonequilibrium population of the spin states 
is only probed by the spin-flip processes. They contribute for eV > B as 

Isf = 27re|jLR|VLPR [{eV - B)P^ + {eV + fi)Pj . (8) 

In equilibrium, only the ground state is occupied, P-^ = 1 and P^ = 0, such that 
only the first term in (8) contributes. Hence, the differential conductance G^^ = 2Gq/3 
remains constant above threshold. In the nonequilibrium situation, the occupation 
probabilities are obtained from the master equation (4) in the stationary state. 




STreljLRpPLPR 



-{eV-B) 
eV -B 




The solution is P^ = 1 — P^ = 1 — 2{eV^B) ■ ^ consequence, now both terms in (8) 
contribute, leading to the total conductance (above threshold) 

G = ho(l + T^^^^] . (10) 



3 V (eV + B)^ 

In the limit V — oo, both P-|-, P^ — )■ 1/2 and the conductance approaches the equilibrium 
value Go with a powerlaw on voltage scale B (although the probability distribution 
remains highly non-equilibrium). 

For the S* = 2 spin of the iron atom with its more complicated spin Hamiltonian, the 
same mechanism as in the simpler spin-1/2 model gives rise to the enhanced conductance 
in the nonequilibrium situation. While in equilibrium only the ground state is occupied, 
Po = 1, leading to steps in the differential conductance, in the nonequilibrium case 
we obtain bias dependent occupations by solving the master equation (4), that lead 
to an overshooting. Again, the conductance decreases above threshold to approach its 
equilibrium value Go in the limit of infinite bias voltage. From our analysis it is clear 
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Figure 2. Fano factor F and dF/dV as a function of bias voltage for different values 
of the asymmetry parameter a. For eV/ksT the Fano factor diverges due to 
thermal noise. Parameters as in figure 1. 



that the nonmonotonic differential conductance is due an increase of transport enabled 
by the population of excited states above threshold but close to the step. It is not a 
signature of the excited spin states carrying less current than the ground state, i.e., a 
decrease of the conductance, as has been speculated in [19]. This conclusion can be 
experimentally checked by measuring the Fano factor, i.e., the ratio between current 
noise and average current, F = S/{eI), as we now explain. 

The Fano factor F is shown in figure 2 for different values of the asymmetry 
parameter a not taking into account the relaxation term (12), see below. For a — )■ 1, we 
find F = 1, i.e., Poissonian behavior, as expected for transport through a normal tunnel 
barrier. When a nonequilibrium population of the atom spin states becomes important 
(a < 1 and bias voltage exceeding the inelastic threshold), the Fano factor becomes 
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super-Poissonian, reaches a maximum and then slowly drops towards the Poissonian 
limit for large bias (The transition |0) — ?■ |1) hardly gives rise to super-Poissonian current 
noise as it is a very weak excitation, cf. figure 1). The latter behavior is an indicator 
that the nonmonotonic conductance is not due to a smaller current contribution from 
the excited states. If this was the case, we would expect a random telegraph signal with 
super-Poissonian Fano factor for — > oo. 

The mechanism leading to the super-Poissonian noise for bias voltages above the 
inelastic thresholds can most easily understood by considering the spin-1/2 model again. 
In this case, the Fano factor above threshold is 

^ . , 2i?^ 3ieV-Bf + 8ieV-B)B 

{eV + BY 3(eV - Bf + 9(eV - B)B + 2B^' ^ ' 

Spin-conserving tunnelling processes are stochastically independent of each other and 
of the spin-fiip transitions. They obey Poissonian statistics and can be ignored for the 
following discussion. Once the inelastic transport channel is open, spin-fiip transitions 
set in. They lead to an alternating sequence of the spin being in the ground and the 
excited state. In the limit of large bias voltage, the rates for the spin-fiip transitions t~^i 
and J,— )-t become equal, and the transport statistics becomes Poissonian. For voltages 
just above threshold, eV > B, however, the two spin-fiip rates differ from each other. 
As a consequence, we obtain an alternating sequence of a longer and a shorter waiting 
time, i.e., effectively there is a tendency of two electrons to bunch together, which yields 
super-Poissonian current noise. 

While our theory predicts an overshooting of the differential conductance at all 
conductance steps, in the experiment of [12] this feature is absent for the steps associated 
with the transition between the ground state |0) and the first excited state |1) whenever 
these steps are pronounced as is the case for a magnetic field along the x axis. This 
indicates that some relaxation mechanism reduces the occupation of We note that 
the transition matrix element of Sz between the ground state and the first excited state is 
large compared to matrix elements of and Sy as well as compared to matrix elements 
of Sz between the ground state and any other excited state. This observation is not 
very sensitive to the direction and the strength of the applied magnetic field. Therefore, 
we make the ad-hoc assumption that there is an additional spin relaxation channel that 
couples to the ^-component of the local spin only. 

We add to our master equation (4) the following phenomenological, spin-dependent 
relaxation rates 

r 

for m 7^ m', where 6(x) is the step function and r is the time scale for relaxation. The 
energy dependence in (12) is not crucial for our conclusions. We therefore choose the 
simplest possible ansatz that allows relaxation only into states with lower energy. In 
contrast, the spin matrix elements are crucial as they suppress the nonequilibrium effects 
for the transition between ground and first excited state while leaving them unaffected 
for almost all other transitions. 
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Figure 3. Differential conductance taking into account a spin-dependent relaxation 
mechanism of the form (12). The relaxation time is given in units of tq with 
= 2TTDS^\jLR\^pi_PR. Parameters are = 3T, T = 0.5K, a = 0. The 
corresponding eigenenergies and eigenstates are summarized in table 1. 



In figure 3 we plot the differential conductance in the presence of a magnetic field 
= 3T along the x-direction. The first transition is more pronounced than in figure 1 
where Bz = 7T along the z-direction. In the limit r — )■ oo we recover the situation 
discussed above where an overshooting effect can be observed for each conductance step. 
By choosing a finite value for the relaxation time comparable to the cotunnelling rates 
exciting the system we can, however, eliminate the overshooting at the first transition 
while leaving the remaining part of the conductance curve practically unaffected. In the 
limit r = 0, we recover the equilibrium value for the conductance at each step. These 
results are not sensitive to the choice of the size and direction of the magnetic field. 

The anisotropic relaxation can also explain the absence of conductance steps due 
to transitions between excited states in the experiment. Such features should be present 
for a small magnetic field applied in the z direction as in this case the first excited state 
gets populated significantly at the first conductance step, the excitation energies satisfy 
A21 < A20 such that the transition |1) — |2) occurs before the onset of the transition 
|0) — )■ 1 2) and furthermore the transition matrix elements (2|S'o|l) do not all vanish. 
However, as the relaxation prevents the system from populating the first excited state, 
these additional conductance features vanish together with the overshooting at the first 
step. 
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4. Conclusions 

We investigated the nonequilibrium effects in transport through a single iron atom. 
With our model, we were able to explain the nonmonotonic features of the differential 
conductance observed experimentally in [12]. Furthermore, we noted a striking absence 
of this nonmonotonicity at certain conductance steps which can be explained by the 
presence of an anisotropic spin relaxation channel. The anisotropy [35] points at the 
importance of spin-orbit coupling in this process, in addition to the splitting of the 
multiplet. In addition, we predicted the occurrence of super-Poissonian current noise 
as a consequence of the nonequilibrium spin occupations probabilities. In conclusion, 
for a full understanding of inelastic tunnelling spectroscopy, it is crucial to account 
for nonequilibrium populations of the atom states established by the competition of 
transport and anisotropic relaxation. 
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